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ABSTRACT 

The collapse of a uniformaly rotating, supermassive star (SMS) to a supermassive black hole (SMBH) 
has been followed recently by means of hydrodynamic simulations in full general relativity. The initial 
SMS of arbitrary mass M in these simulations rotates uniformly at the mass-shedding limit and is 
marginally unstable to radial collapse. The final black hole mass and spin have been determined to be 
Mh/M ss 0.9 and Jh/M% « 0.75. The remaining mass goes into a disk of mass M^sk/M w 0.1. Here 
we show that these black hole and disk parameters can be calculated analytically from the initial stellar 
density and angular momentum distribution. The analytic calculation thereby corroborates and provides 
a simple physical explanation for the computational discovery that SMS collapse inevitably terminates 
in the simultaneous formation of a SMBH and a rather substantial ambient disk. This disk arises even 
though the total spin of the progenitor star, J/M 2 = 0.97, is safely below the Kerr limit. The calculation 
performed here applies to any marginally unstable n — 3 polytrope uniformly rotating at the break-up 
speed, independent of stellar mass or the source of internal pressure. It illustrates how the black hole and 
disk parameters can be determined for the collapse of other types of stars with different initial density 
and rotation profiles. 

Subject headings: black hole physics - relativity - hydrodynamics - stars: rotation 



1. INTRODUCTION 

Recent observations provide increasingly strong evi- 
dence that supermassive black holes (SMBHs) of mass 
~ 10 6 — lO lo Af0 exist and that they are the central en- 
gines that power active galactic nuclei (AGNs) and quasars 
(Rees 1998, 2001). The dynamical formation of SMBHs, 
as well as the inspiral, collision and merger of binary 
SMBHs, are promising sources of long-wavelength gravita- 
tional waves for the proposed Laser Interferometer Space 
Antenna (LISA) (Thorne 1995; Schutz 2001). However, 
the scenario by which SMBHs form is still very uncertain 
(see, e.g., Rees 1984, for an overview). One promising- 
route is the collapse of a supermassive star (SMS). Once 
they form out of primordial gas, sufficiently massive stars 
will evolve in a quasistationary manner via radiative cool- 
ing, slowly contracting until reaching the point of onset 
of relativistic radial instability. At this point, such stars 
undergo catastrophic collapse on a dynamical timescale, 
leading to the formation of a SMBH (Bisnovatyi-Kogan, 
Zel'dovich & Novikov 1966; Zel'dovich & Novikov 1971; 
Shapiro & Teukolsky 1983). 

Because most objects formed in nature have some an- 
gular momentum, rotation is likely to play a significant 
role in the quasistationary evolution, as well as the final 
collapse of a SMS. The slow contraction of even a slowly 
rotating SMS will likely spin it up to the mass-shedding 
limit, because such stars are so centrally condensed. At 
the mass-shedding limit, matter on the equator moves in 
a Keplerian orbit about the star, supported against gravity 
by centrifugal force and not by an outward pressure gra- 



dient. Baumgarte & Shapiro (1999) recently performed 
a detailed numerical analysis of the structure and stabil- 
ity of a rapidly rotating SMS in equilibrium. Assuming 
the viscous or magnetic braking timescale for angular mo- 
mentum transfer is shorter than the evolution timescale 
of a typical SMS (Zel'dovich & Novikov 1971; Shapiro 
2000), the star will settle into rigid rotation and evolve 
to the mass-shedding limit following cooling and contrac- 
tion. They found that all stars at the onset of quasi-radial 
collapse have an equatorial radius R » 640GA/ / c 2 and a 
nondimensional spin parameter a/M = cJ/GM 2 m 0.97. 
Here J, M, c, and G are spin, gravitational mass, light 
velocity and gravitational constant. (Hereafter we adopt 
gravitational units and set c = G = 1). Because of the 
large value of a/M, it was uncertain whether the rotating 
SMS would collapse directly to a black hole and/or form 
a disk. Because of the growth of T"/|W| oc 1/R during col- 
lapse (where T is the rotational kinetic energy and W is 
the gravitational potential energy) , it was unclear whether 
the collapse would trigger the growth of bars or other non- 
axisymmetric instabilities. We therefore investigated the 
collapse of a rotating SMS with a 3 + 1 post-Newtonian 
(PN) hydrodynamic simulation (Saijo et al. 2002). We 
found that the collapse proceeds nearly homologously and 
axisymmetrically and inferred that a SMBH is likely to be 
formed, with some gas remaining outside the hole. Shi- 
bata & Shapiro (2002) then performed simulations of the 
collapse using an axisymmetric code in full general rela- 
tivity. We showed conclusively that a black hole forms at 
the center and determined that the mass of the hole con- 
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tains about 90% of the total mass of the system and has 
a spin parameter J/M 2 ~ 0.75. The remaining gas forms 
a rotating disk about the nascent hole. 

Guided by our numerical simulations,we show here that 
the final black hole mass and spin can be calculated an- 
alytically (up to quadrature) from the stellar density and 
angular momentum profiles in the progenitor SMS. The 
calculation provides a simple physical explanation for the 
important finding that a SMBH formed from the collapse 
of a rotating SMS is always born with a "ready-made" 
disk that can provide fuel for accretion to power an en- 
ergy source. 

2. BASIC ASSUMPTIONS 

Our analysis relies on several explicit assumptions, all 
of which we expect to hold to high accuracy in SMSs. In 
particular, we assume that the SMS is 

1. dominated by thermal radiation pressure; 

2. fully convective; 

3. uniformly rotating; 

4. governed by nearly Newtonian gravitation; 

5. described by the Roche model in the outer 
envelope. 

For large masses, the ratio between radiation pressure, 
P r , and gas pressure, P g , satisfies 
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(see, e.g., Shapiro & Tcukolsky 1983, eqs. (17.2.8) 
and (17.3.5)); here the coefficient has been evaluated for 
a composition of pure ionized hydrogen. For SMSs with 
M <; 10 6 M Q , we can therefore neglect the pressure con- 
tributions of the plasma in determining the equilibrium 
profile, even though the plasma can be important for de- 
termining the stability of the star when it is rotating suffi- 
ciently slowly (Zel'dovich & Novikov 1971; and Shapiro & 
Teukolsky 1983). A simple proof that SMSs are convective 
in this limit is given in Loeb & Rasio (1994). This result 
implies that the photon entropy per baryon, 

4 aT 3 

S ^3— (2) 
is constant throughout the star, and so therefore is /? w 
8(s r /k)~ 1 . Here a is the radiation density constant, n is 
the baryon density, and k is Boltzmann's constant. As a 
consequence, the equation of state of a SMS is that of an 
n = 3 polytrope: 
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P = Kp 4 /\ K = 



= const, (3) 



where m is the atomic mass unit and /x is the mean molec- 
ular weight (cf. Clayton 1983, eq. 2-289; note that Clayton 
adopts a different definition of (3, which is related to ours 
by /Sciayton = /?/ (1 + In fact, our results depend only 
on the assumption that the star is an n = 3 polytrope, 
and not on the origin of the internal pressure. 

The third assumption, that the star is uniformly rotat- 
ing, is probably the most uncertain of our assumptions. 



Nevertheless, it has been argued that convection and mag- 
netic fields provide an effective turbulent viscosity which 
dampens differential rotation and brings the star into uni- 
form rotation (Bisnovatyi-Kogan, Zel'dovich & Novikov 
1967; Wagoner 1969; Shapiro 2000; for an alternative, see 
New and Shapiro 2001). Uniformly rotating configura- 
tions were adopted as initial configurations in Shibata and 
Shapiro (2002), whose numerical results we are attempting 
to corroborate. 

We assume that gravitational fields in the progenitor 
star are sufficiently weak so that we can apply Newto- 
nian gravity in analyzing its envelope. SMSs at the on- 
set of collapse have equatorial radii R/M w 640 (Baum- 
garte & Shapiro 1999), so that this assumption certainly 
holds. Relativistic corrections are crucial for the stabil- 
ity of SMSs, but can be neglected in the analysis of the 
equilibrium state. 

Finally, the Roche approximation provides a very ac- 
curate description of the envelope of a rotating stellar 
model with a soft equation of state, as in the case of 
an n — 3 polytrope (for numerical demonstrations, see, 
eg, Papaloizou & Whelan 1973 and Baumgarte & Shapiro 
1999). Since our analysis is based on this approximation, 
we will briefly review it together with some of its predic- 
tions in the following section. 

We assume that the collapse proceeds 

1 . axisymmetrically ; 

2. adiabatically. 

Our previous 3D simulation in PN gravitation (Saijo et 
al. 2002) allowed for the growth of nonaxisymmetric per- 
turbations, but none were observed. Hence the assumption 
of axisymmetry, which was the basis of the fully relativis- 
tic simulation of Shibata & Shapiro (2002), appears to be 
justified. Cooling by photon and neutrino emission can be 
ignored on a dynamical timescalc (Linke et al 2001), hence 
the collapse is essentially adiabatic. The collapse will be 
accompanied by a burst of gravitational radiation which 
may be detectable by LISA. However, the fractional loss of 
mass-energy is sufficiently small that it can be ignored in 
determining the black hole mass (Shibata & Shapiro 2002); 
angular momentum is strictly conserved, as gravitational 
waves carry off no angular momentum in axisymmetry. 

3. REVIEW OF THE ROCHE MODEL 

Stars with soft equations of state are extremely centrally 
condensed: they have an extended, low density envelope, 
while the bulk of the mass is concentrated in the core. 
For an n = 3 polytrope, for example, the ratio between 
central density to average density is p c /p = 54.2. The 
gravitational force in the envelope is therefore dominated 
by the massive core, and it is thus legitimate to neglect the 
self-gravity of the envelope. In the equation of hydrostatic 
equilibrium, 

— = -V($ + $ c ), (4) 
P 

this neglect amounts to approximating the Newtonian po- 
tential $ by 

M 

* = "— ( 5 ) 
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In (4) wc introduce the centrifugal potential <& c , which, 
for constant angular velocity f2 about the z-axis, can be 
written 



$ c = -htf (x 2 + y 2 ) = -\tt 2 r 2 sin 2 9. 
Integrating eq. (4) yields the Bernoulli integral 

h + $ + $ c = H, 

where H is a constant of integration and 



h = 



dP . 1 P 
— = n+1 - 
P P 



(6) 



(7) 



(8) 



is the enthalpy per unit mass. Evaluating eq. (7) at the 
pole yields 



H = - — , 

Rn 



(9) 



since h = on the surface of the star and <I> C = along 
the axis of rotation. In the following we assume that the 
polar radius R p of a rotating star is always the same as 
in the nonrotating case. This assumption has been shown 
numerically to be very accurate (e.g., Papaloizou & Whe- 
lan 1973). The mass of the star is hardly changed from its 
value in spherical equilibrium, 

M = ^(^'^)l)^ 3/2 ' ( 10 ) 

where the Lanc-Emdcn factor for an n = 3 polytrope is 
given by £ 2 |0'(£i)| = 2.01824 (Shapiro & Teukolsky 1983). 
A rotating star reaches mass shedding when the equator 
orbits with the Kepler frequency. Using eqs. (6) and (7), 
it is easy to show that at this point the ratio between 
equatorial and polar radius is 



L P / shedd 

The corresponding maximum orbital velocity is 



a 



shedd 



2\ 3/2 /MA 



(11) 
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(Zel'dovich & Novikov 1971; Shapiro & Teukolsky 1983). 

Inserting eqs. (3), (5), (6), (8), (9) and (12) into eq. (7) 
yields the density throughout the extended envelope, 

(£|0'(£i)|) 2 M (R p 1 , 4 r 2 ^ 2a \ 3 (13) 



P = 
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R* \ r 



4 r 2 , 

1 + -sin 2 6> 
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The scale of the density is set by the value of the polar 
radius at the point of onset of radial collapse, 

R p /M = 427 (14) 

(Baumgarte & Shapiro 1999). The stellar surface is the 
boundary along which p = 0, and is thus defined by the 
curve r{9) given by 

4r3 -- " 1 = 0. (15) 



27 Rf" 2 6 R~ p 



The solution to this cubic equation is given by 
r(9) _ 3 sin (0/3) 



Rn 



sin 9 



(16) 



c!*>q 




(R. Cooper, private communication). 



Fig. 1. — A SMS rotating uniformly at the mass-shedding 
limit. The shaded region shows the portion of the Roche enve- 
lope ra > roisco where the matter has specific angular momentum 
j > ilSCO- Consequently, this gas will escape capture by the black 
hole that forms at the center following collapse of the interior. 

4. CALCULATING THE BLACK HOLE MASS AND SPIN 

Consider the implosion of matter from the envelope onto 
the central black hole formed from the collapse of the cen- 
trally condensed interior region. If the specific angular mo- 
mentum of the imploding envelope matter is below JiscO) 
the specific angular momentum of a particle at the in- 
nermost stable circular orbit (ISCO) about the hole, the 
matter will be captured. If the angular momentum of the 
infalling matter exceeds jisco, the matter will escape cap- 
ture and continue to orbit outside the hole, forming a disk. 
This capture criterion is well-supported by the numerical 
simulation of Shibata & Shapiro (2002), where the cap- 
ture of fluid in noncircular orbits with j > jisco was 
demonstrated to be negligible. This criterion suggests a 
simple iterative scheme for calculating the final mass and 
spin of the hole and disk from the initial stellar density 
and angular momentum profile. First, guess the mass and 
spin of the hole, Mh and J^. For our initial guess, we 
shall take a black hole that has consumed all the mass 
and angular momentum of the star, so that Mh/M = 1 
and Jh/Ml = 0.97. Next, use the initial stellar profile 
to "correct" this guess by calculating the escaping mass 
and angular momentum of the outermost envelope with 
specific angular momentum exceeding jisco (see figure 1). 
We note that the value jisco depends on Mh and Jh- We 
then "correct" the black hole mass and spin by deducting 
the values of the escaping mass and angular momentum of 
the envelope material from the guessed values of Mh and 
Jh- We recompute jisco f° r the "corrected" black hole 
mass and spin, and repeat the calculation of the escaping 
envelope mass and angular momentum until convergence 
is achieved. The calculation described above relies on the 
theorem that for an axisymmetric dynamical system, the 
specific angular momentum spectrum, i.e., the integrated 
baryon rest-mass of all fluid elements with specific angular 
momentum j less than a given value (e.g., jisco) is strictly 
conserved in the absence of viscosity (Stark & Piran 1987). 
Any viscosity, if present, is expected to be unimportant on 
dynamical timescales, as required by the theorem. 

For a Kerr black hole of mass Mh and spin parameter 
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a = Jh/M h , the value of jisco is given by 



Jisco 



VM h r ms (rl s - 2a^M h r ms + a 2 ) 



• ms 

(rl s - 3M h r ms + 2a^/M h l~) 1 / 2 
where r ms is the ISCO given by 

r ms = M h {3 + Z 2 - [(3 - Zi)(3 + Z 1 + 2Z 2 )]^ 2 }, 
where 



a 



and 
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(17) 



(18) 



1/3" 



(19) 



(20) 



(see, e.g., Shapiro & Teukolsky 1983). Clearly, the infalling 
gas corotates with the black hole. 

The mass of the escaping matter in the envelope with 
3 > Jisco is given by 



AM 



-II 



2-Kwdwp 



(21) 



where the density is given by eq. (13) and the quadrature is 
performed over all cylindrical shells in the star with cylin- 
drical radii w > wisco = (jisco/^) 1 ^ 2 (see fig 1). Here 
we set f2 = fi s hedd given by eq. (12). Defining m = vj/R p 
and z = z/Rp gives the nondimensional integral 

AM/M = (£i|#'(£i)|) 2 

3 



// 



zudzudz 



1 



[w 2 + Z 2 )!/ 2 



(22) 

Similarly, the angular momentum carried off by the escap- 
ing matter in the envelope is given by 



AJ 



2irzudTXjp Sim 2 



(23) 



or 



AJ/M 2 = (£|0'(£i)|) S 



3/2 



1/2 



// 



m 3 dmdz 



M 

1 _£_ 2 

{w 2 +z 2 y/ 2 27 W 



1 3 



(24) 

where R p /M is given by eq. (14) and where once again the 
integral is performed over all cylindrical shells in the star 
with w > wisco = (j'isco/^) 1/2 

The mass and angular momentum of the black hole can 
then be determined from eqs. (22) and (24) according to 

M h /M = 1 - AM/M (25) 

and 

( J/M 2 - AJ/M 2 ) 



Jh/M 2 = 



(26) 



(1 - AM/M)' 

Once the iteration of eqs. (17) - (20) with cqs. (22), (25) 
and (26) converges, the mass of the ambient disk can be 
found from 

M disk /M = AM/M. (27) 

Convergence to better than 1% is achieved after only 
four iterations. We find that vuisco/Rp = 0.43. We ob- 
tain a black hole mass Mh/M — 0.87 and black hole spin 



J h /Ml = 0.71, and a disk mass M disk /M = 0.13. These 
values are consistent with those inferred from the numeri- 
cal simulation of Shibata & Shapiro (2002): M h /M w 0.9, 
Jh/Ml w 0.75, and M disk /M ss 0.1. The differences are 
well within the numerical accuracy of the dynamical inte- 
grations in the simulation. 



5. DISCUSSION 

We have considered a SMBH formed by the collapse 
of a radially unstable, uniformly rotating SMS of mass 
M spinning at the mass-shedding limit. We calculated 
that the black hole has a mass Mh/M 0.87 and a spin 
Jh/Ml w 0.71. Most significant, we found that at birth 
such a black hole will be automatically embedded in a sub- 
stantial disk of mass M^isk/M w 0.13. A disk forms even 
though the total angular momentum of the star at the 
onset of collapse is J/M 2 w 0.97, below the Kerr limit. 

We have also calculated the black hole and disk param- 
eters for unstable SMSs which, at the onset of collapse, 
are spinning more slowly than the angular frequency at 
break-up (i.e. mass-shedding). Wc find that the parame- 
ters are very insensitive to the ratio a = f2/Sl s h c dd, where 
^shedd is given by eq. (12), over the range 0.1 < a < 1. 
The reason is that for small f2, a rotating SMS destabi- 
lizes at large R p /M (Rp/M w 427/a 2 ) with a value of 
J/M 2 w 0.97 nearly independent of a. The correspond- 
ing value of tnisco /Rp ~ 0.43 is nearly unchanged, and 
as a result, the final black hole and disk parameters do 
not vary significantly as a decreases. (For sufficiently low 
spin rates, the effects of thermal gas pressure compete 
with rotation to determine the onset of stability and must 
be included). Nevertheless, we believe, on evolutionary 
grounds, that it is most likely that a SMS will be spinning 
near the break-up value at the onset of collapse (Baum- 
garte & Shapiro 1999). 

The nascent disk is fairly massive. Such a massive, ro- 
tating disk may be dynamically unstable to nonaxisym- 
mctric instabilities. If the disk deforms, it could emit long 
wavelength, quasi-periodic gravitational waves with a fre- 
quency of about 10- 3 (M/10 6 M Q )- 1 Hz, which might be 
detectable by LISA. 

At birth, the disk will be hot and thick, radiation 
dominated, and will extend from radius ~ 640M at the 
outer edge in the equatorial plane all the way down to 
the horizon. Subsequently, the matter in the disk will 
likely cool by both photon and neutrino emission and flat- 
ten. Under the influence of viscosity, the gas will diffuse 
outwards in the outermost region, transporting angular 
momentum away from the inner region. At the same 
time, gas in the inner region will accrete onto the cen- 
tral black hole, increasing its mass and powering a source 
of luminous energy. The lifetime of the disk tdisk may 
be estimated crudely by speculating that the accretion 
eventually will settle into a steady state, generating pho- 
tons at the Eddington luminosity at ~ 10% efficiency: 
L = L Ed d = 1.3 x 10 38 (M/j/Mq) and M = 10L/c 2 . Ac- 
cordingly, we have t d isk = M disk /M, or 



^disk — 



M disk c 2 

10L E dd 



5 x 10 6 yrs. 



(28) 
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The above estimate for idisk is independent of M. How- 
ever, it may be off by a considerable factor: a lower lu- 
minosity will increase the disk lifetime, while a lower ef- 
ficiency (e.g, an efficiency characterizing an ADAF disk; 
Narayan & Yi 1994) or an enhancement in the total lumi- 
nosity via the addition of neutrino cooling, could decrease 
the lifetime. Moreover, the disk at black hole birth may 
comprise far more material than the gas which escapes 
capture during the final implosion. Prior to undergoing 
collapse, a rotating SMS will likely evolve along a mass- 
shedding sequence as it cools and contracts, provided the 
stellar viscosity and/or magnetic braking proves sufficient 
to enforce uniform rotation (Baumgarte & Shapiro 1999). 
In this case, the secular cooling and contraction phase will 
be accompanied by mass loss, and the gas which is ejected 
then might already form an extensive, ambient disk prior 
to collapse. Finally we note the possibility that such a 
massive disk around the hole may suffer a global "run- 
away" instability that may cause a more rapid destruc- 
tion of the disk than implied by steady-state accretion 
(Abramowicz et al. 1983; Nishida et al. 1996; Nishida 
& Eriguchi 1996). These are important issues which are 
ripe for further analysis. 

The calculation performed here applies to any 
marginally unstable n — 3 polytrope rotating at the mass- 
shedding limit. The results are independent of stellar mass 



or the source of internal pressure. The method should be 
applicable for estimating the black hole and disk param- 
eters for the collapse of other types of stars with differ- 
ent equations of state, initial density and rotation profiles. 
For example, a preliminary calculation indicates that for 
a marginally unstable star uniformly rotating at break-up 
speed, the mass fraction that avoids collapse and forms a 
disk falls sharply as the polytropic index of the star drops 
below n = 3. The reason is that such stars are more 
compact at the onset of collapse, with smaller values of 
Rp/M and, hence, less extended envelopes. This may pro- 
vide a qualitative explanation for the results of previous 
simulations in full general relativity for collapse of rotat- 
ing neutron stars modeled as n = 1 polytropes (Shibata, 
Baumgarte & Shapiro, 2000; Shibata, 2000). We hope to 
provide further examples in the future. 
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of Illinois at Urbana-Champaign and Monbu-Kagaku-sho 
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